Growth and properties of strained VO^ thin films with controlled stoichiometry 
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We have succeeded in growing epitaxial films of rocksalt VO x on MgO(OOf) substrates. The 
oxygen content as a function of oxygen flux was determined using 18 02-RBS and the vanadium 
valence using XAS. The upper and lower stoichiometry limits found are similar to the ones known 
for bulk material (0.8<x<1.3). From the RHEED oscillation period a large number of vacancies 
for both vanadium and oxygen were deduced, i.e. wl6% for stoichiometric VO. These numbers 
are, surprisingly, very similar to those for bulk material and consequently quite strain-insensitive. 
XAS measurements reveal that the vacancies give rise to strong low symmetry ligand fields to be 
present. The electrical conductivity of the films is much lower than the conductivity of bulk samples 
which we attribute to a decrease in the direct overlap between t2 g orbitals in the coherently strained 
layers. The temperature dependence of the conductivity is consistent with a variable range hopping 
mechanism. 



I. INTRODUCTION 

One of the major challenges in transition-metal oxide 
thin film research is to grow the material in single crystal 
form with controlled oxygen stoichiometry. This is not a 
trivial task especially when the transition metal ion has 
multiple valences and the variation of the oxygen con- 
tent leads to a rich and complex phase diagram. Impor- 
tant examples for such materials include the binary oxide 
systems like the vanadium and titanium oxides 1-5 . The 
most common and convenient growth methods in thin 
film research like pulsed laser deposition and sputtering 
techniques are not very suitable here, because of the high 
substrate temperatures and consequently also high oxy- 
gen pressures usually used, giving less opportunity to fine 
tune the oxygen stoichiometry. These techniques are ob- 
viously the method of choice for ternary and quarternary 
oxides, in which one can vary the material properties by 
tuning one of the constituent cation concentrations while 
keeping that of the oxygen fixed. For the binary oxides, 
however, one has to resort to true molecular beam epi- 
taxy (MBE) techniques, in which one has to carefully 
balance the oxygen flux with respect to the metal evap- 
oration rates in order to obtain good control and tuning 
of the stoichiometry. The question here is whether one 
can obtain thin films with physical properties that are 
systematic, controllable and tunable. 

Here we report on our research project to grow single 
crystal thin films of vanadium monoxide (VO x ), as well 
as on the study of the physical properties and the elec- 
tronic structure of these films. Bulk vanadium monoxide 
has many intriguing properties that are closely related to 



the issue of stoichiometry. While at first sight the crys- 
tal seems to have the very simple fee rocksalt structure, 
it actually can have a very wide range of varying oxy- 
gen concentrations: values of x = 0.8 — 1.3 have been 
reported for bulk VO x . Even more remarkable is that 
there is a large number of both cation and anion vacan- 
cies, even for a; = 1, in which case the concentrations 
are as high as 15% 6 . Stoichiometric VO always remains 
disordered, whereas ordering of the vacancies was only 
reported for x between 1.2 and 1.3 7 . A systematic inves- 
tigation of the physical properties of polycrystalline VO x 
has been carried out by Banus et al. 6 in the early 70's. 
VO x shows an interesting transition from a metallic to 
a semiconducting type of behavior at x close to 1, and 
the magnetic susceptibility can be described by a Curie- 
Weiss law, suggesting that the V 3d electrons partially 
localize 3 . Several models have been proposed to explain 
these properties, in which the role of electron correlation 
effects, band formation and vacancies are discussed 8-11 . 

Only very little work on films of VO^ has been done 
up to know. Metal supported vanadium oxides were 
investigated by several authors in connection with cat- 
alytic properties. Thin V oxide layers deposited on 
Cu(100) 12 ' 13 and Ni(110) 14 were found to consist of a 
V203-like oxide phase after room temperature oxidation. 
VO^-like islands growth was reported on Rh(lll) 15 . 
Surnev et al. found a VO/VO2 - like phase at low cov- 
erage on Pd(lll) using scanning tunnelling microscopy 
(STM) 16 ' 17 . Thin films of VO x have also been grown 
on Ti02 substrates, because of the enhanced catalytic 
activity shown by titania-supported vanadium oxides 18 . 
Sambi et al. 19 ~ 21 studied locally ordered VO^ films with a 
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thickness up to 4 ML grown on TiO2(110) by evaporating 
vanadium metal at room temperature and subsequent an- 
nealing in vacuum for a short period of time, in order to 
promote the uptake of oxygen from the substrate. From 
a structural point of view, local order of their films is 
supported by X-ray photoelectron diffraction (XPD). 

In the examples mentioned above, the vanadium oxi- 
dation state was established on the basis of the position 
and shape of the V 2p 3 / 2 photoemission line. However, 
the uncertainty of stoichiometry of vanadia films remains, 
since the binding energies (BE) of V 2p 3 / 2 for different ox- 
idation states are very close, a comparison between data 
obtained in Refs. 12 and 19 showing some discrepancies. 
Moreover, a determination of the physical properties of 
VOa, thin films is missing in these studies. Preparation 
is often found to be difficult since VOa; easily oxidizes 
to V2O3. These problems can be avoided by growing 
epitaxial films of VO^ on appropriate substrates with a 
careful optimization of the evaporation conditions. 

This paper describes how we have been able to suc- 
cessfully grow coherent epitaxial single crystalline VOa; 
thin films by stabilizing them on MgO(lOO) substrates. 
The paper is organized as follows. The details of 
the experimental conditions and equipment used arc 
given in the experimental section (section II). The re- 
sults are described in several subsections. We first 
present (section IIIA) the growth process as monitored 
structurally in- situ by reflection high-energy electron 
diffraction (RHEED) and low-energy electron diffrac- 
tion (LEED) and compositionally by x-ray photoelectron 
spectroscopy (XPS). We then (section IIIB) determines 
the stoichiometry x of the various VO x films by using 
Rutherford backscattering spectrometry (RBS). Further 
structural details of the epitaxy (section IIIC) is revealed 
by ex-situ x-ray diffraction (XRD). The important is- 
sue of vacancy concentrations is addressed using a com- 
bination of the RBS and RHEED intensity oscillations 
(section HID). Basic data about the electronic structure, 
such as the valence and local crystal fields of the ions, 
(section HIE) are measured using soft-x-ray absorption 
spectroscopy (XAS). We end the results section by pre- 
senting the temperature dependent transport properties 
of the films (section IIIF). Finally we discuss the results 
by comparing the properties of the VO x thin films with 
those of the bulk material, and with the predictions made 
by existing theoretical models . 



II. EXPERIMENTAL 

The experiments were performed in an ultra-high vac- 
uum (UHV) MBE system, with a base pressure below 
IxlCT 10 mbar, equipped with RHEED, LEED and XPS 
facilities. RBS and XRD were used to characterize the 
thin films grown in the MBE system ex-situ. MgO(lOO) 
was chosen as the substrate because it has the same 
cubic rocksalt structure as VO x , with a mismatch of 
about 3%. MgO blocks were cleaved ex-situ and immedi- 



ately introduced into the UHV chamber where the sub- 
strates were cleaned by annealing at 650° C for 2 hours in 
lxlO -6 mbar oxygen to remove hydrocarbon contamina- 
tions. Clean and well-ordered surfaces, as determined by 
XPS, RHEED and LEED were obtained. Vanadium (V) 
(purity 99.99%) was evaporated using an electron-beam 
evaporator (Omicron EFM3). The deposition rate of V, 
measured by moving a quartz crystal at the position of 
the substrate, was set to 0.7 — 0.8 A/min. The vanadium 
was simultaneously oxidized by a beam of 18 02 (99.99%). 
The oxygen was supplied through a leak valve into a 
small buffer volume 22 connected to the vacuum cham- 
ber through a 35 cm long and 1 cm wide stainless steel 
pipe, ending at a distance of 10 cm from the substrate. 
The pressure in the buffer volume was measured with 
a Baratron capacitance manometer. The variation in 
buffer pressure during deposition was less than 1%. VO x 
films with different stoichiometries were grown by vary- 
ing the 18 02 buffer volume pressure, while keeping the V 
flux constant. During evaporation the background pres- 
sure was 2xl0 -8 mbar or lower. The V and 18 C>2 beams 
can be abruptly and simultaneously stopped. After clos- 
ing the 18 2 valve the background pressure dropped to 
below 10 -9 mbar within seconds, insuring a well-defined 
oxygen exposure of the sample. All the samples were 
grown while keeping the substrate at room temperature 
(RT). 

RHEED was used to monitor the evolution of the films. 
Thicknesses were calibrated by monitoring the RHEED 
intensity oscillations during deposition. Oscillations in 
the RHEED specular beam intensity, where each oscil- 
lation corresponds to the formation of one new atomic 
monolayer (ML), allows for control of the film thickness. 
Film thicknesses were also determined ex-situ from X- 
ray specular reflectivity (XRR) measurements 23 ' 24 and a 
good agreement with RHEED intensity oscillations was 
obtained. 

XPS measurements were performed using non- 
monocromatic Al Ka radiation (hv = 486.6 eV) and 
the total energy resolution of the electron analyzer in 
combination with the photon source is about 1 cV. To 
avoid ex-situ post-oxidation, a thin MgO cap layer (20A) 
was grown for protection by Mg deposition in 18 02 at- 
mosphere of lxlO -8 mbar at room temperature. Con- 
sequently all the ex-situ measurements, i.e. RBS, XRD, 
XAS and electrical measurements were done on capped 
VOa; samples. 

To determine the oxygen content, RBS measurements 
have been performed. For calibration purposes a V2O3 
epitaxial thin film was also grown on an Al 2 03(0001) 
substrate, with the oxygen partial pressure and the 
substrate temperature set at 10~ 6 mbar and 550°C 
respectively 25-27 . All the measurements were done in a 
so-called 'random' orientation with respect to the crys- 
talline axes, using a He + beam with 1.5 MeV energy. 
The backscattering yields were detected at an angle of 
165 degrees. 

The XRD measurements were done using a Philips 
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FIG. 1: RHEED patterns recorded at an electron energy of 
15 kV with the beam incident along [100] direction, (a) clean 
MgO(lOO) substrate.(b) 20A VO X (100) thin film grown on 
MgO(lOO) substrate using 0.18 mTorr oxygen buffer pressure 
while keeping the substrate at room temperature. The picture 
is inverted in order to observe better the Kikuchi lines. 

MRD diffractometer, equipped with a hybrid mir- 
ror/monochromator for Cu K Q radiation, a 4-circle go- 
niometer and a programmable slit in front of the detector. 

The X-ray absorption spectroscopy (XAS) measure- 
ments were performed at the 11 A Dragon beamline 28,29 
at the Synchrotron Radiation Research Center Taiwan, 
using the total electron yield mode. The light has a de- 
gree of linear polarization of 98% and an energy resolu- 
tion of 0.16 eV for photon energies between 500 and 550 
cV. The spectra are normalized to the photon flux mea- 
sured using a gold mesh. The angle of the incident light 
is normal to the VO x films. 

Resistivity measurements were performed in a com- 
mercial Quantum design PPMS system. The electrical 
measurements were done in the constant voltage mode 
and the current was measured along the [100] direction 
of the films by two and four point geometry. Electrical 
contacts consisting of 10 nm of Cr metal were evaporated 
ex-situ on the MgO substrates prior to introduction into 
the MBE system, in order to measure the resistance of 
as-grown VO x films. 



III. RESULTS 
A. Growth of VO^ thin films 

For a wide range of incident oxygen fluxes, epitax- 
ial growth of VOz thin films on MgO(100) substrates 
was concluded from in situ RHEED and LEED anal- 
ysis. Our MBE system geometry allows for monitor- 
ing changes in surface structure during film growth by 
RHEED. This low angle electron diffraction technique 
is particularly suited to thin film growth as it is highly 
sensitive to surface morphology. A RHEED pattern of 
a clean MgO(100) surface and 20A thick VO^ thin film 
deposited on MgO(100) is shown in Fig.l. 

The beam was incident along the [100] direction. The 
basic RHEED pattern did not change during growth, sug- 
gesting that the film continues to grow as a rocksalt phase 
on top of the underlying substrate. Sharp streaks and the 
presence of Kikuchi lines indicate that the surface is still 
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FIG. 2: RHEED intensity oscillations of the specularly re- 
flected electron beam observed during deposition of VO^ on 
MgO(100). The electron beam was incident along the [100] 
direction, with a primary energy of 15 kV. 



smooth after deposition of 10 monolayers (ML) of VO^. 

In most cases oscillations in the intensity of the specu- 
larly reflected beam were observed during growth. These 
oscillations are characteristic for a layer-by-layer growth 
mode, each oscillation corresponding to the formation of 
one atomic layer. Note that RHEED oscillations were 
observed at room temperature, which is considered to be 
a very low temperature for an oxide system. In Fig. 2 an 
example of these growth oscillations is shown. 

In Fig. 3 a sequence of RHEED patterns which were 
taken after deposition of about 50 monolayers of VO x 
using different oxygen fluxes is presented. 

The oxygen content x is indicated on each RHEED pat- 
tern. The corresponding oxygen flux used can be found 
in the caption of Fig. 3. In paragraph III.B we will explain 
how the x values were determined using RBS spectrom- 
etry. Varying the oxygen buffer pressure between 0.12 
and 0.20 mTorr (0.82 < x < 0.97) epitaxial growth is 
obtained and sharp streaks can still be observed after 
deposition of 50 ML VCv The distance between VO^ 
streaks is the same as between the MgO streaks. This is 
consistent with a coherent growth as determined ex-situ 
by XRD, which will be discussed in section III.C. The 
rocksalt structure of VO x thin films was also confirmed 
by LEED. The LEED pattern displays the same square 
symmetry and periodicity as a MgO(100) surface. 

Between 0.22 and 0.30 mTorr (1.1 < x < 1.31), 
MgO RHEED streaks fade away quickly after starting 
the growth and 3D-transmission spots appear, suggesting 
considerable roughening of the surface. The absence of a 
LEED pattern for the x = 1.28 and x = 1.31 cases con- 
firms that the surface morphology is becoming increas- 
ingly disordered. Nevertheless, the presence of RHEED 
oscillations shows that the growth is still layer-by-layer 
like. 

For an oxygen buffer pressure of less than 0.12 mTorr 




FIG. 3: Sequence of RHEED patterns observed after deposi- 
tion of 50 ML VO x on MgO(lOO) at different oxygen fluxes 
while keeping the substrate at room temperature. The oxygen 
buffer pressure was varied between 0.10 and 0.30 mTorr: (a) 
0.30 mTorr, (b) 0.27 mTorr, (c) 0.24 mTorr (d) 0.22 mTorr, 
(e) 0.20 mTorr, (f) 0.18 mTorr, (g) 0.16 mTorr , (h) 0.14 
mTorr, (i) 0.12 mTorr, (j) 0.10 mTorr. The electron beam 
was incident along [100] direction, with a primary energy of 
15 kV. 



(x < 0.8) powder rings and extra spots were observed 
in RHEED. Vanadium atoms were not completely oxi- 
dized in the later case as deduced from XPS core level 
intensities. 

Attempts to make VO x films at higher substrate tem- 
perature were unsuccessful. The RHEED pattern quickly 
disappeared or was getting very diffuse from the begin- 
ning. 

The \-2p and O-ls core lines measured for three dif- 
ferent films i.e vanadium metal on MgO(100), VO x on 
MgO(100) and V 2 3 on Al 2 O 3 (0001) are shown in Fig.4. 
The binding energies were corrected for the charging ef- 
fect by assuming a constant binding energy of the Ols 
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FIG. 4: V 2p and O Is XPS spectra of (a) V 2 3 on A1 2 3 
(0001), (b) VOz on MgO(OOl) and (c) vanadium metal on 
MgO(OOl). 



peak at 531 eV. All spectra were corrected in a standard 
manner for the satellites due to the Ka^ay components of 
the incident X-ray and an integral background was sub- 
tracted afterwards. The binding energies of the V 2p 3 / 2 
peak found are given in the figure. The binding energy 
corresponding to the V 2p 3 / 2 in the VO x film was found 
between the ones corresponding to V metal and V2O3, re- 
spectively (see Fig.4) . The vanadium oxidation state can 
be in principle determined from the position and shape 
of the V 2p3/2 feature. Nevertheless, due to the fact that 
vanadium with different oxidation states has rather sim- 
ilar values for the binding energies 30 we did not perform 
a quantitative analysis of the XPS data. In order to de- 
termine precisely the stoichiometry we have adopted a 
different route, which will be presented in the following 
section. 



B. Stoichiometry determination 

The oxygen content of the films was determined from 
RBS measurements. As was mentioned already in the 
experimental part, 18 O2 instead of 16 O2 was employed 
for film growth to distinguish between the oxygen of the 
film and substrate. Using a 1.5 MeV He + beam, a good 
mass separation between 16 O from the substrate and 18 O 
from the film can be obtained. This can be observed in 
Fig. 5 showing a RBS spectrum of a V2O3 layer grown 
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FIG. 5: RBS spectra of 1.5 MeV He + ions scattered from a 
V2O3 film epitaxial grown on an Ai2O3(0001) substrate. The 



O peak can be clearly distinguished. 



FIG. 7: The oxygen content 1 as a function of the buffer 
volume pressure as determined from RBS data for two series 
of VOi samples. 



VO capped 



"'0 



VO x uncapped 



200 220 240 260 280 

Back scattering energy (keV) 



300 



FIG. 6: RBS spectra of 1.5 MeV He ions scattered from 
VO x film deposited on MgO(100). (a) with a thin (2nm) 
MgO caplayer, (b) without caplayer. An extra shoulder at 
the 16 O position is observed due to postoxidation. 



epitaxially on an Al 2 O3(0001) substrate, which was used 
for calibration. A nice hexagonal LEED pattern charac- 
teristic for the corundum structure was observed, proving 
the high quality and long range order of the V2O3 film. 
In Fig. 6 part of the spectra for a capped and uncapped 
VOz layer on MgO(OOl) is shown. In both spectra a well- 
separated 18 0-peak is visible. For the uncapped layer, an 
additional shoulder at the 16 position can be observed 
due to the fact that the layer is post-oxidized in air. The 
absence of this shoulder in the spectrum of the capped 
layer proves that capped VO x films are really protected 
from post-oxidation. 



With the following formula we have calculated the ratio 
between the numbers of V and O atoms in a VO x film, 
using V2O3 as a reference sample. 



3 [ a°] vo x 
lA v >v 2 o 3 



(1) 



The areas Ao and Ay in the RBS spectra were de- 
termined by fitting the V and 18 O peaks from VO x and 
V2O3 samples to a gaussian function, in the case of the 
18 0-peak after subtraction of the linear background due 
to Mg. The 18 0-peak area was corrected for the contri- 
bution of the caplayer, using the thickness of the VO^ 
films and the Mg ls O caplayer, as determined from the 
RHEED oscillations. The resulting s-values are plotted 
in Fig. 7 as a function of the oxygen buffer volume pres- 
sure for two series of samples. One can observe that the x 
values are nicely reproducible. The error in the determi- 
nation is about 3%. The upper and lower stoichiometry 
limit found are similar to the ones known for bulk mate- 
rial. For the sample grown at 0.10 mTorr the vanadium 
metal is not completely oxidized, as deduced from the 
XPS spectra. For x bigger than 1.3 a mixture of two 
phases was found. 



C. Structure of VOj thin films 

From RHEED and LEED patterns it was already con- 
cluded that VO x films grow epitaxially on MgO(100). To 
analyze the epitaxy and crystal structure of the films in 
a more quantitative way, we also performed an ex-situ 
XRD analysis. 

The measurements were done on samples capped with 
a thin (20A) epitaxial MgO film. 9-29 scans show only 
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FIG. 8: 8-20 X-ray diffraction measurement of lOnm thick 
VO x thin film prepared with x = 0.94. Only the (002) and 
(004) diffraction peaks characteristic of the rocksalt structure 
can be observed. The reflections from the film are broadened 
due to the finite thickness of the film. Subsidiary thickness 
fringes indicate a well defined composition and thickness of 
the film. 



weak diffraction peaks close to the (002) and (004) peaks 
of the MgO substrate, as expected from films having a 
rocksalt structure. The reflections are broadened due 
to the finite thickness of the film. There are no peaks 
corresponding to other phases. Moreover, for samples 
made with the oxygen buffer pressure varying between 
0.12 and 0.22 mTorr a number of subsidiary thickness 
fringes were observed, suggesting a well denned compo- 
sition and thickness of the film (see Fig. 8). 

To determine whether the growth is fully coherent or 
partially relaxed, the most convenient way is to mea- 
sure the intensity profile around a non-specular reflection 
common to the substrate and the overlayer. This can be 
done by performing a set of 20-ui scans at different uj, 20 
and uj being the detector and sample orientations with 
respect to the beam direction. These scans can be easily 
mapped into reciprocal space. An example of such a map 
for the region in reciprocal space around the (113) reflec- 
tion is shown in Fig. 9. The horizontal and vertical axes 
are k vectors parallel (k par ) and perpendicular (k per ) to 
the surface plane, respectively. The intensity scale in the 
figure is logarithmic. The feature corresponding to the 
VO x film can be clearly distinguished. The elongated 
shape of the MgO reflection perpendicular to the radial 
direction is due to mosaic spread. The peaks of MgO 
and VO x are at the same k par value, which proves that 
the film is fully coherent. Consequently, VO^ thin films 
experience a compressive strain in the perpendicular di- 
rection, which is induced by the lattice mismatch. 

The lattice constant normal to the substrate surface 
can be determined from — 20 scans around MgO (002) 
or MgO (004) peaks using the Bragg law. 20 values were 
determined by fitting the peaks to a gaussian function. 
The perpendicular lattice constant is plotted in Fig. 10 as 
a function of the oxygen content. 

We found that the perpendicular lattice constant and 
consequently also the average lattice constant is decreas- 
ing with increasing oxygen content. This is in contrast 
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FIG. 9: XRD map of reciprocal space around the non-specular 
(113) reflection of a lOnm thick VOj, thin film with x = 0.94. 
The logarithm of the diffracted intensity as a function of the 
in-plane k» and out-of-plane k±_ reciprocal lattice vectors is 
plotted. The x- and y- axis are in units of 3?, with A = 
0.15015 nm. 
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FIG. 10: Perpendicular lattice constant of VOz films with 
different oxygen content as a function of x as determined from 
XRD. The perpendicular lattice constant is decreasing with 
increasing x. 



to the behavior reported by Banus 6 for bulk VO x . 

The layer thicknesses determined from X-ray specular 
reflectivity (XRR) measurements were in good agreement 
with the values obtained from the RHEED intensity os- 
cillation period. 



D. Vacancy concentration 

As was already mentioned in the Introduction, a large 
number of both vanadium and oxygen vacancies is char- 
acteristic for bulk VO x . Using a combination of RHEED 
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FIG. 11: Vanadium and oxygen vacancy concentrations in 
VOi as determined from the period of RHEED oscillations. 
The number of vanadium vacancies is increasing with x, while 
the number of oxygen vacancies is decreasing. 



and RBS resuls, an estimation of the number of vacancies 
present in VO^ thin films can be made. RHEED oscilla- 
tion periods were determined for all the VO x samples as 
well as for a calibration V2O3 thin film on AkiC^OOOl). 
All samples were grown using the same value for the vana- 
dium flux. By comparing the time needed to grow one 
monolayer of VO x and V2O3 and assuming that V2O3 
is stoichiometric and free of vacancies, the fraction of 
vacant vanadium sites (Vy) in VO^ can be directly cal- 
culated. The oxygen vacancy concentration (Vo) can be 
obtained because x is already known from RBS, using 
the following expression: 



l-Vo 
l-V v 



(2) 



As illustrated in Fig.ll, the vanadium vacancy con- 
centration is increasing with increasing oxygen content 
x, while the oxygen vacancy concentration is decreasing. 
The total number of vacancies (Vy+Vo) is also decreas- 
ing with x. Close to the stoichiometric value we found 
about 16 % vacant sites of both vanadium and oxygen. 



E. Valence of the V ions from XAS 

XAS measurements were done at the O K- edge in 
order to verify the valence of the V and O ions. The 
data are depicted in Fig. 12 for samples with the oxygen 
content x (as determined from RBS) between 0.8 and 
1.3. For comparison, we have also included the spec- 
trum of a V2O3 film in top part of the figure. The O K- 
edge absorption spectra correspond to the dipole-allowed 



transitions from the O Is to the O 2p shell, which is par- 
tially empty due to the hybridization with V 3c? conduc- 
tion band states. The spectral structures that can be 
observed in the 528 — 534 eV photon energy range are 
dictated by these V 3d states 31-34 . For photon energies 
higher than 535 eV, the oxygen of the MgO cap layer also 
start to contribute to the XAS signal 34 . 

The distinct peak observed at photon energies around 
532 eV, can be assigned to transitions into the empty 
V 3d-e g ^ band 33 . Structures at higher photon energies 
can be ascribed to transitions to the higher lying V 4sp 
related bands. Important is to note that for the most 
oxygen deficient VO x samples, i.e. x « 1, the lowest 
spectral structure is given by the 532 eV peak, indicating 
that the lower lying V 3d-t2 g ^ band is completely occu- 
pied by three electrons, and that thus the V valence is 2+ 
or less. For samples with higher oxygen content, i.e. for 
x > 1, a clear low energy peak appears at about 530 eV. 
This strongly suggests that holes are introduced in the V 
3d-t2 g ^ band, meaning that the V valence is higher than 
2+. This assignment is supported by the V2O3 spec- 
trum, in which transitions to both the and the e g ^ 
bands are possible because only two electrons occupy the 
three-fold degenerate orbital in this V 3+ 3c? 2 system. 

These XAS measurements show that the cross-over 
from less than 2+ to more than 2+ V valencies occurs for 
an RBS x value of about 0.94 — 0.97. This value is not 
very far from 1.00, and can be taken as an indication for 
a good agreement between the XAS and RBS methods. 

We have also carried out XAS experiments at the V 
£23 (2p — > 3c?) edges. The results are shown in Fig. 13, in 
which we have also included the spectra for V2O3 and 
Cr203 for reference purposes. In going from x = 0.8 to 
1 = 1.3 we can observe a gradual change in the spectra. 
Distinct and sharp structures start to develop for x > 1. 
The similarity of these structures with those of V2O3 
is striking, and in fact indicates that V ions in a 3+ 
valence state are present for x > 1, consistent with the 
observations at the O K edge mentioned above. 

Important is the observation that the spectra are quite 
broad for x around 1. In the simplest approximation, one 
would expect to see a spectrum with the typical atomic 
multiplet structure of a 2p 6 3c? 3 — » 2p 5 3c? 4 transition of 
a V 2+ ion in a Oh symmetry, as shown in the bottom 
curve of Fig. 13. This curve has been calculated using 
standard parameters (10D q ~ 2 eV) often applied for 
the analysis of transition metal oxide soft-x-ray absorp- 
tion spectra 35-37 and it shows quite distinct features with 
several peaks and valleys. Clearly, the experimental V 
L23 spectra for x ~ 1 do not have these distinct features. 

Also a comparison with the experimental spectrum of 
Cr2C>3 (see Fig. 13), which is also a 3d 3 ion in approxi- 
mately Or local symmetry, leads to the conclusion that 
the V L23 spectra are anomalously broad. We take this 
observation as an indication that the vanadium in VO is 
not all in the local 3c? 3 — Oh symmetry, but that instead 
an appreciable amount of vanadium ions are experienc- 
ing strong local ligand fields of low symmetry associated 
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FIG. 12: O Is X-ray absorption spectra of VO^ samples w 
the oxygen content varying between 0.8 and 1.3. 



with the presence of large amount of vacancies. Th( 
low symmetry ligand fields must be at least several hx 
dreds of meV strong, in order to wash out completely t 
multiplet structure of a 3c£ 3 -ion in Oh symmetry. lm 
implication of such fields will be discussed in Section IV. 



F. Transport properties 

Electrical measurements have been performed to study 
the electronic structure of VO^ films. Resistivity data 
have been reported previously by Banus 6 for polycrys- 
tallinc samples. Bulk material exhibits a semiconducting 
behavior for x > 1, with an activation energy rising to 
about 40 meV for x = 1.3. For x < 1 VO x behaves 
like TiO x , with an almost temperature and composition 
independent resistivity of about 3xl0 -3 ilcm. Banus et 
al. reported a transition from a semimetallic to semicon- 
ducting type behavior in VO x at x = 1.05. 

The variation of the logarithm of the electrical resis- 
tivity of VOz thin films with different oxygen content is 
shown as a function of 1/T in Fig. 14. The resistivity is 
increasing with decreasing T for all but the x = 0.82 sam- 
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FIG. 13: V L23 X-ray absorption spectra of VOx samples with 
the oxygen content varying between 0.8 and 1.3. 



pie, suggesting a semiconductor-like behavior in the en- 
tire temperature range studied. However, it is evident in 
Fig. 14 that logp does not vary linearly with 1/T. Fig. 15 
shows the composition dependence of the activation ener- 
gies calculated from the logp-vs-l/T plots close to room 
temperature. Note, that the apparent activation energy 
decreases with decreasing temperature. 

Such behavior is similar to that described by Banus 6 , 
but we did observe a few important differences. Com- 
pared to the corresponding resistivity results for bulk 
samples we found a much higher absolute value of the 
resistivity. As can be observed from Fig. 14 the transi- 
tion from a metallic to a semiconductor type behavior is 
shifted from x = 1, as found in bulk material, to x = 0.8 
for thin films. Moreover, the activation energy calculated 
close to room temperature (Fig. 15) is larger. For exam- 
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FIG. 14: Logarithm of resistivity of VOi films with different 
oxygen content plotted against 1/T. Resistivity is increasing 
with decreasing temperature except for the x = 0.82 sample, 
suggesting a semiconductor-like behavior. 
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FIG. 15: Activation energy E a versus x calculated for the 
temperature interval 150 < T < 300 K. 



pie, for bulk stoichiometric VO, Banus et al. 6 obtained an 
apparent activation energy of only 5 meV. For thin films 
we found a much higher value of « 110 meV. A large 
increase in the activation energy can be observed going 
from x = 0.94 to x = 0.97. Although the oxygen content 
is changing only with 3%, the activation energy is be- 
coming almost twice as large. For x « 1.3 the activation 
energy rises to 180 meV. 

The T-dependence of the resistivity is better described 
by variable range hopping behavior 38 given by: 



FIG. 16: Logarithm of normalized resistivity of VO x thin 
films with different oxygen content plotted as a function of 

T -l/4 



TABLE I: Parameters extracted from T 1//4 dependence of 
resistivity. 

xinVO, T (K) |N(E F )(cm- 3 eV- 1 ) 



0.82 
0.84 
0.89 
0.94 
0.97 
1.1 
1.28 



2.6 
2xl0 3 
5xl0 4 
12xl0 4 
5xl0 6 
lxlO 7 
2xl0 7 



2xl0 26 
2xl0 23 
8xl0 21 
3xl0 21 
lxlO 20 
4xl0 19 
2xl0 19 



/T i 
P ~ exp(y)4 



(3) 



This is evident from Fig. 16 showing the normalized re- 
sistivity versus T -1 / 4 . This result is not surprising know- 
ing that a high number of vacancies is characteristic for 
this material. A random field caused by vacant sites can 
produce localization of the electrons near the Fermi level 
just like in amorphous materials. Its occurrence in VOa; 
seems to be strong evidence that the vacancies play a 
major role, causing large changes in its electronic prop- 
erties. The parameter To can be extracted from the slope 
of the curves in Fig. 16. T is related to the density of lo- 
calized states at the Fermi level N(E^) 39 by the following 
expression: 



To 



24a J 



nN(E F )k B 



(4) 



where cv _1 is the decay length of the wave function as- 
sociated with the charge carriers and ks is the Boltz- 
mann constant. Taking a reasonable value for cv _1 as 5A, 
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N(Eir) can be calculated from the values of To. These 
parameters are tabulated in Table 1 together with the 
corresponding x values. Clearly, one can observe that for 
x = 0.82 sample the density of localized states at the 
Fermi level is high and starts to decrease with increasing 
the oxygen content. 



IV. DISCUSSION 

The transport properties of VO^ can not be under- 
stood on the basis of a simple band picture, because the 
material should have been metallic. Shortly after the 
publication of Banus et al. , two papers by Goodenough 8 
and Mott 9 appeared, in which explanations were given in 
terms of electron correlation effects and vacancies. Both 
authors assume a finite Hubbard U, which is large enough 
to split the t 2g band. Because V 2+ supplies 3 electrons 
to the t 2g band, the valence band of VO is given by a 
completely filled lower Hubbard band, separated by a 
gap from the conduction band which is the empty upper 
Hubbard band. 

According to Mott 9 , the random field produced by the 
vacancies causes Anderson localization 40 in the overlap 
region of the two Hubbard bands. This explains why 
the conductivity in VO x is of the variable range hop- 
ping type at low temperatures. In Goodenough's model 8 
the effect of the vacancies on the electronic structure is 
discussed much more explicitly. He assumes a high de- 
gree of trapping of electron and holes near anion and 
cation vacancies, respectevely. In this way the Fermi- 
level stays near the minimum in the DOS of VO^, ex- 
plaining its semiconductor-like behavior. The transition 
to a more itinerant behavior for x < 1 is attributed to a 
tail in the trapped electron distribution. A second effect 
of the trapping of charges near vacancies is that the loss 
of Madelung energy is limited. 

The presence of a large number of vacancies in VO x re- 
mains one of the most puzzling characteristics. The only 
two other binary oxides showing this behavior are the 
rocksalts TiOa; and NbO, with NbO having about 25% 
vacancies 41 . Apparently the creation of vacancies stabi- 
lizes the crystal by reducing the Gibbs free energy. Ac- 
cording to Goodenough 8 the formation of vacancies leads 
to a reduction in the lattice constant, thereby broadening 
the t 2g bands. The resulting stabilization of the occu- 
pied valence band states is counteracted by a reduction 
in Madelung energy, but this energy loss is assumed to 
be minimized by a localization of the charge compensat- 
ing electrons and holes near the oxygen and vanadium 
vacancies respectively. He argues that the application of 
hydrostatic pressure should reduce the vacancy concen- 
tration considerably and he refers to experiments giving 
indirect evidence in favor of this effect 42 . 

The VO x films we have grown epitaxially on MgO sub- 
strates are thin enough to be coherent with the substrate 
and are under tensile stress. In fact, the in-plane lat- 
tice constant is expanded (~3%), while the out-of-plane 



lattice constant is reduced due to the Poisson effect. Ap- 
plying Goodenough's vacancy induced lattice contraction 
arguments, we would expect that the formation of vacan- 
cies will not be less favored in our thin film, because the 
lattice constant of the film is fixed by the substrate, thus 
preventing the energy gain to occur that otherwise is as- 
sociated with the lattice contraction and broadening of 
the tig bands. In contrast to this expectation, however, 
we found that both the cation and anion concentration 
are not much affected. For stoichiometric VO film, the 
values are even slightly higher (of the order of 1 or 2%) 
than for bulk material. 

We agree with Goodenough, nevertheless, that the 
mechanism for vacancy formation must be searched for 
in terms of energy arguments. To our knowledge, ab- 
initio total energy calculations have not yet been car- 
ried out to study the stability of vacancy formation in 
this material. We may speculate about the mechanism 
by reviewing Goodenough's arguments 8 . Goodenough 
starts with the picture that electrons and holes are be- 
ing trapped near the cation and anion vacancies, as to 
avoid part of the Madelung energy loss. The t 2g orbitals 
of cations neighboring a cation vacancy are then desta- 
bilized by a reduced bonding to neighboring t 2g orbitals 
and a stronger tt bonding to the p orbitals next to the va- 
cancy, thereby raising their energy above the Fermi level. 
On the other hand, near an oxygen vacancy the bonding 
t 2g orbitals are stabilized, whereas the e g and s orbitals 
are less destabilized by the absence of oxygen p a orbitals 
at the vacant site. These effects contribute to the stabi- 
lization energy of vacancies, but not enough. Instead of 
looking for a possible further energy gain in the t 2g band 
due to lattice contraction induced band broadening, we 
now discuss what low symmetry ligand fields can do to 
lower the energy of the occupied t 2g , motivated by the 
fact that the XAS measurements reveal that those ligand 
fields are very strong in VO^, i.e. with a strength of the 
order of several hundreds meV. 

We will first discuss the case of oxygen vacancies. 
Neighboring divalent 3d metal ions are surrounded by five 
oxygen ligands in a square pyramidal arrangement with 
C4 V symmetry. The t 2g and e g levels will be split into 
three non- and one doubly degenerate levels, respectively. 
With two electrons being trapped at or near the vacancy, 
we have one extreme possibility in that both electrons are 
trapped at the vacant site, forming an F-center, and the 
other possibility in that the excess electrons are trapped 
at neighboring ti-metal sites. Assuming that at least a 
fraction of the excess electrons resides at the (i-metal 
ions, the average number of d electrons will be larger 
than three for V-ions next to a vacancy. In both cases 
some ligand field stabilization energy is gained, owing to 
the electrons occupying the lower e g level in VO x <i. 

Considering now the case of d metal vacancies, we may 
speculate that two excess holes are trapped near these 
vacancies in VO^i. One extreme possibility is that the 
holes are residing at oxygen neighbors. Following a sug- 
gestion of Elfimov et al. 43 , this would give a stable config- 
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uration in which the holes are in a triplet state localized 
at the oxygen coordination polyhedron around the va- 
cancy. However, in contrast to the case of CaO discussed 
by these authors, in a transition metal compound there is 
also the possibility that the two holes are located at next 
nearest neighbor transition metal ions. The presence of 
the cation vacancy leads to a local symmetry lowering 
(C2t,) at these cation sites. Also in this case a ligand field 
stabilization is anticipated, because the average number 
of d electrons on ci-metal ions next to a vacancy will be 
less than three for V. 

Although, as was argued above, the decrease of direct 
overlap between t^g orbitals in the strained films does not 
have a major effect on the concentration of vacancies, it 
will change the band structure in the valence band re- 
gion and consequently change the electrical properties. 
Expanding the lattice by epitaxial growth on a substrate 
with a larger lattice constant like MgO will decrease the 
bandwidth and therefore increase the gap. This explains 
the lower conductivity and higher room temperature ac- 
tivation energies found in our films. In the models pro- 
posed by Goodenough 8 and Mott 9 , the on-site Coulomb 
interaction U is assumed to open up a small energy gap 
in the itinerant ti g band. It is assumed that overlapping 
tails of localized states, which are related to the disorder 
associated with the large number of defects in the sys- 
tem, are present at the edges of these Hubbard bands. 
In this way Mott 9 explains why in bulk samples the con- 
ductivity vs temperature is semiconductor-like at high 
temperature, but is better described by a variable range 
hopping mechanism at low temperatures. The conduc- 
tivity results for our films are consistent with this ex- 
planation. The main difference is that the resistivity and 
apparent activation energies of the films are much higher, 
and that the transition from semiconducting to metallic 
behavior is shifted to lower x values, i.e. from x w 1 in 
bulk samples to x ~ 0.8 in films. 



V. CONCLUSIONS 

We have successfully grown epitaxial VO^ films on 
MgO(lOO) substrates. Up to at least 120A the growth 
is coherent and layer- by-layer-like. 18 02-RBS was intro- 
duced as a convenient method to determine accurately 
the stoichiometry of these ultrathin layers. Once the sto- 
ichiomctry is known, the vacancy concentration of both 
vanadium and oxygen can be calculated from the time 
to grow one monolayer as determined from RHEED. The 
numbers turn out to be very similar to those for the bulk 
material. This implies that the formation of a high con- 
centration of vacancies may not be directly related to 
an increased ti g band width as a result of a vacancy 
induced lattice contraction. Instead we suggest that a 
detailed study is required to calculate the possible sta- 
bilization due to additional ligand field splittings at the 
low symmetry metal sites near the vacancies, as observed 
from XAS. Nevertheless, the decrease in direct overlap 
between ti g orbitals and the concomitant increase of the 
size of the pseudo-gap between the lower and upper ti g 
Hubbard bands is held responsible for the much larger 
electrical resistivity in strained VO x films. 
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